Direct observation of charge inversion by multivalent ions as a universal electrostatic 

phenomenon 
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We have directly observed reversal of the polarity of charged surfaces in water upon the addition of 
tri- and quadrivalent ions using atomic force microscopy. The bulk concentration of multivalent ions 
at which charge inversion reversibly occurs depends only very weakly on the chemical composition, 
surface structure, size and lipophilicity of the ions, but is dominated by their valence. These results 
support the theoretical proposal that spatial correlations between ions are the driving mechanism 
behind charge inversion. 



Understanding screening due to mobile ions in liquid 
is a key theme of such diverse fields as polymer physics, 
nanofluidics, colloid science and molecular biophysics. 
Screening by multivalent ions in particular results in sev- 
eral counter-intuitive phenomena, for example attraction 
between like-charged macromolecules such as DNA l] 
and actin filaments '2| . Similarly, the electrophoretic mo- 
bility of charged colloids reverses sign upon introducing 
a sufficient concentration of multivalent ions in solution 
[3|, 1^ , a phenomenon known as charge inversion. 

The conventional paradigm for describing screening in 
liquid divides the screening ions into two components: 
(1) the so-called Stern layer, consisting of ions confined to 
the surface, and (2) a diffuse component described by the 
Poisson-Boltzmann (PB) equation that decays exponen- 
tially with distance far from the charged object. Charge 
inversion can be accounted for by introducing a "chem- 
ical" binding constant that reduces the free energy of 
multivalent ions situated in the Stern layer, reflecting an 
assumed specific interaction between these ions and the 
surface being screened. This chemical binding constant is 
expected to depend on properties of the ions such as their 
size, chemical composition, surface structure, lipophilic- 
ity and valence. While this approach has been successful 
in describing experimental data 0j IM IE Hi i^ usually 
provides little insight into the underlying binding mech- 
anism and lacks significant predictive power. 

A universal mechanism for charge inversion based pre- 
dominantly on electrostatic interactions has been pro- 
posed Q- It was noted that the predicted chemical po- 
tential of the Stern layer can be significantly lowered if 
spatial correlations between discrete ions are accounted 
for. At room temperature, the loss of entropy entailed 
by the formation of a highly correlated ionic system is 
substantial. For multivalent counterions and sufficiently 
high surface charge densities, however, this is more than 
compensated by the corresponding gain in electrostatic 
energy, leading to charge inversion 9]. To date, these 
theories have remained untested by experiments. 

Here we present direct measurements of charge inver- 
sion and its dependence on the properties of the screen- 
ing ions. Using an atomic force microscope (AFM), we 



measured the force between two oppositely charged sur- 
faces. This approach circumvents the main limitations of 
previous measurements, namely, reliance on modelling of 
hydrodynamic effects 0, 3 ^tnd the need to disentangle 
phenomena at two similarly-charged surfaces |^ J^ . We 
observe that in the presence of a sufficiently high concen- 
tration of tri- and quadrivalent ions, the force reversibly 
changes sign. The bulk concentration at which charge 
inversion occurs, the so-called charge-inversion concen- 
tration cg, depends almost exclusively on the valence of 
the ions, consistent with the universal predictions of the 
ion-correlation theories. 

Positively charged amine-terminated surfaces were pre- 
pared under argon atmosphere by immersing silicon 
wafers with 200-500 nm thermally-grown oxide in a 0.1% 
solution of 1-trichlorosilyl-ll-cyanoundecane (Gelest) in 
toluene for 30 minutes, then in a 20 % solution of Red 
Al (Sigma-Aldrich) in toluene for 5 hours. Negatively 
charged surfaces were prepared by gluing 10 /im diam- 
eter silica spheres (G. Kisker Gbr) with epoxy resin to 
AFM cantilevers (ThermoMicroscope Microlevers) using 
a method similar to that of Ducker et al [lOj, as illus- 
trated in Fig. ^a). 

Using a Digital Instrument NanoScope IV AFM, force 
spectroscopy measurements were performed yielding the 
force F on the silica bead versus the bead-surface sepa- 
ration d [lO| . The spring constant of the cantilevers was 
0.03 N/m, as given by the manufacturer. Care was taken 
to minimize the scattering of light from the surface so as 
to eliminate interference effects. 

At separations d greater than the Debye length A of 
the solution, the force decays exponentially with d: 



i^ = i^oexp(-rf/A), d>X. 



(1) 



The parameter Fq is proportional to the so-called renor- 
malized surface charge densities of both the silica bead 
and the aminc-terminated surface, a^ and a* respec- 
tively. The values of crj ^ are related by the PB equation 
to the net surface charge densities ab and CTs (including 
both the bare surface charge and the charge in the Stern 
layer). At low net surface charge densities \(Jb,s\ < Cmax, 
the renormalized charge densities are simply equal to the 



net charge densities: cr^g = <7b,s- Here amax = 4 kTe/eX, 
where k is the Boltzmann constant, T is the temperature, 
e is the dielectric constant of water and — e is the elec- 
tron charge. At higher net charge densities, a^ ^ saturates 
at (Tniax- Because we use oppositely charged surfaces and 
Z:l electrolytes, where Z is the valence of the multivalent 
ions, only one of the surfaces is affected by the introduc- 
tion of multivalent ions. The other surface thus plays 
the role of a constant probe. Near charge inversion, Fq is 
thus proportional to the net surface charge density of the 
surface being screened by multivalent ions, ab or as, and 
the sign of the force yields unambiguously the polarity of 
this net surface charge. 

For d ^ A, the PB equation predicts a more compli- 
cated form than Eq. (^. In addition, van der Waals 
forces, regulation of the surface charge and depletion 
forces can become important. We therefore concentrate 
our analysis on the regime d > X. 

Three positive trivalent ions. Lanthanum La'^"'", ruthe- 
nium(III) hexammine [Ru(NH3)6]^^ and cobalt(III) 
sepulchrate [CoCi2H3oN8]'^+ were investigated. La'^+ 
is an elemental metal ion with a first hydration shell 
consisting of 8-9 water molecules (radius r of the com- 
plex 398 pm milllllEi). [Ru(NH3)6]3+ contains a 
Ru(III) core with six NH3 groups around it (r — 364 pm 
m El Ei)- [CoCi2H3oN8]3+ is a caged cobalt com- 
plex with CH ^ gr oups exposed to the water environment 
(r = 445 pm 15]), making it less hydrophillic than the 
other two. 

Figure ^ shows the measured force-distance relation 
F(d) as a function of multivalent ion concentration c 
for the multivalent salts LaCl3 (b), CoCi2H3oNgCl3 (c) 
and Ru(NH3)6Cl3 (d). A force measurement with only 
supporting electrolyte (LaCl3: [la], C0C12H30N8CI3 and 
Ru(NH3)6Cl3: [13) was first performed (black squares), 
showing an attractive interaction between the surfaces. 
Solutions with increasing concentrations of multivalent 
ions in the monovalent supporting electrolyte were then 
pumped through the AFM fluid cell of 50 /zl volume at 
a rate 0.15-0.2 ml/min for at least 5 minutes per solu- 
tion. This allowed the surface to equilibrate with the 
electrolyte and insured that c was not influenced when 
large numbers of ions screened the surface. Consecutive 
measurements of F(d) at multivalent ion concentrations 
c = 10 /xM, 100 /zM and 1 mM are shown in Fig.[l] At the 
end of the experiment, the measurement with c = 10 fiM 
was repeated (red open circles). The C0C12H30N8CI3 
and Ru(NH3)6Cl3 measurements were carried out con- 
secutively using the same silica bead. 

We interpret these observations as follows. The posi- 
tive multivalent ions adsorb on the negative silica bead, 
reducing <7b and thus the magnitude of the force. Near 
1 mM, the screening charge in the Stern layer overcom- 
pensates for the bare surface charge; at becomes positive 
and the force becomes repulsive. The last measurement 
with c = 10 /iM, which shows a recovery to the force mea- 
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FIG. 1: (color) (a) Optical microscope images of the side 
(left) and top (right) of a cantilever with a silica sphere. Force 
versus separation measurements in different concentrations 
of (b) LaCla, (c) CoCiaHsoNsClg (c) and (d) Ru(NH3)6Cl3. 
Insets illustrate schematically the attractive (1) and repulsive 
(2) forces between the silica bead and the amine-terminated 
surface. The legend applies to all three graphs. 



sured at the beginning of the experiment, indicates that 
charge inversion reflects reversible equilibrium between 
the surface and the bulk electrolyte. 

To further compare the charge-inversion concentration 
of the same surface with different multivalent ions, the 
force for d > A was fitted to Eq. (^ . Because it is difficult 
to accurately fit the Debye length A when the force is 
very small, its value was fitted for the curve with c — 
and corrected using the standard expression for A for the 
cases c > 0. 

Figure [2Ia) shows the fitted normalized force extrap- 
olated to zero separation, Fjvo(c) = Fo(c)/Fo(0), for the 
[CoCi2H3oN8]3+ and [Ru(NH3)6]3+ data of Fig. [Ti;c,d). 
Similarly, Fig. ^h) shows Fjvo(c) for consecutive mea- 
surements using the same silica bead on La'^"'' (data from 
Fig. mb)) and [Ru(NH3)6]3+ (F(d) curves not shown). 
We estimate the charge- inversion concentration cq by lin- 
early interpolating between the data points immediately 
above and below Fn = on the lin-log scale. In both sets 
of measurements, the observed values of co differ by a fac- 
tor ~ 2. More generally, we find that the charge- inversion 
concentrations of silica for the three chemically different 
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FIG. 2; Normalized force extrapolated to zero separation 
obtained from fits to Eq. (0, versus multivalent ion concen- 
tration c for (a) C0C12H30N8CI3 (squares) and Ru(NH3)6Cl3 
(circles) and for (b) La^+ (squares) and Ru(NH3)6Cl3 (cir- 
cles). In each case the data were obtained consecutively using 
the same silica bead. Lines are guides to the eye. 



trivalent ions La3+, [Ru(NH3)6]3+ and [CoCi2H3oN8]3+ 
differ by at most a factor of 2.1, as summarized in Ta- 
ble m This is comparable to the variation observed be- 
tween measurements for the same ion and pH using dif- 
ferent, nominally identical beads and surfaces. Although 
the charge-inversion concentrations of the three positive 
trivalent ions are similar, there are differences in the ob- 
served F{d) curves. In particular, La'^^ is less effective 
in reducing the absolute force at low concentrations, but 
the magnitude of the force for c ^ cq is largest for La'^+. 

Figure 131 shows measurements where the same amine- 
terminated surface was consecutively charge inverted by 
a molecule in two different charge states, iron(II) hexa- 
cyanide [Fe(CN)6]''^ {r = 443 pm) and iron(III) hcxa- 
cyanide [Fe(CN)6]^" (r = 437 pm) jHElI^J, ensuring 
that essentially the only difference between the two mea- 
surements is the valence of the ions |l8| . Figure |3Ic) 
shows the Fjvo for both ions as a function of the concen- 
tration. The charge-inversion concentrations for the two 
ions differ by a factor ~ 50. 

Measurements using [Fe(CN)6]'^^ and ruthenium(II) 
hexacyanide [Ru(CN)6]'^" (r = 456 pm HI El El E3), 
two ions with nearly identical chemical groups exposed 
to solution and differing only by their core atom, gave 
nearly identical F{d) curves at all concentrations. 

Two divalent ions, calcium Ca^+ and magnesium Mg^+ 
(radii of 388 and 348 pm, respectively ElIllElEi) did 
not show charge inversion at a concentration of 1 mM on 
a silica bead that showed charge inversion at 1 mM La'^+. 
Thus divalent ions, if they can charge invert a silica bead 
at all, do so at higher concentrations than trivalent ions. 
Concentrations higher than ImM were not investigated 
in this study because the Debye length then becomes so 
short that effects such as van der Waals forces mask the 
electrostatic interaction between the surfaces. 

Additional experiments were performed with positively 
charged surfaces made by chemically modifying a sil- 
icon dioxide surface with 3-aminopropyltriethoxysilane 
(APTES) and by adsorbing poly-L-lysine on mica. Key 



0.2 
0.0 
-0.2 
-0.4 
-0.6 



(a) 


^ (b) 


/ 


^■■^^ 


77 

• ■ 


p^^^ 


[Fe(CN),]'- ■ 

- OnM 
. luM 

- 3nM 

- lOuM 


[Fe(CN),]3-- . 
OnM 

. lOnM 
- loonM 

T 500nM 



20 40 60 80 100 120 20 40 60 80 100 120 

d^m) d (nm) 



S 0.0 




10 100 1000 

c(iM) 

FIG. 3: (color) Force versus separation measurements in dif- 
ferent concentrations of (a) K4Fe(CN)6 and (b) K3Fe(CN)6. 
(c) Normalized force at zero separation versus multivalent ion 
concentration c for K4Fe(CN)6 (squares) and K3Fe(CN)6 (cir- 
cles). Lines are guides to the eye. 



results are summarized in Tabled 

In terms of a chemical binding description, our mea- 
surements indicate that the binding constants for La^^, 
[Ru(NH3)6]3+ and [CoCi2H3oN8]3+ on silica differ by 
at most a factor ^ 2, despite the fact that these ions 
have significantly different chemical composition, surface 
structure, size and lipophilicity. The binding constant 
differs at least 10-fold for the same molecule in two dif- 
ferent charge states on amine-terminated surfaces. These 
observations strongly suggest that specific chemical in- 
teractions are not responsible for charge inversion in our 
measurements and that the mechanism for adsorption is 
predominantly electrostatic. 

We compare our results with ion-correlation theories 
using the formalism of Shklovskii '3|, in which the mul- 
tivalent counterions in the Stern layer are assumed to 
form a strongly correlated liquid with short-range corre- 
lations resembling those of a Wigner crystal. This theory 
provides a simple analytical prediction for the charge- 
inversion concentration: 

Co = ((Tbarc/2erZ) cxp{-fi,/kT) exp(-A//fcr). (2) 

Here (Xbaie is the bare surface charge density, A^'^ is the 
standard energy of adsorption of an ion and ^c is the 
chemical potential of the strongly correlated liquid. The 
latter can be approximated by the value for a Wigner 

crystal: fic « ^bare^'^^^- ^"^ ^^^ calculations we use the 
full expression for ^c |9|. In the absence of hydration 
effects and specific chemical interactions, A^° = and 
fic is the sole driving force behind charge inversion. 

Qualitatively this theory with A/i° = is in good 
agreement with our observations. The predictions that 



TABLE I: Summary of measurements in which the same surface was charge inverted by two different ions. 



surface 


probe 


supp. elect. 


ion(l) 


ion(2) 


4" (PM) 


4'' (mM) 


(high) 1 (low) 


APTES 


sihca bead 


m 


[Fe(CN)6]*- 


[Fe(CN)6]3- 


13 


170 


13 


chlorosilane 


sihca bead 


m 


[Fe(CN)6]^- 


[Fe(CN)6]3- 


4 


200 


50 


APTES 


sihca bead 


\M 


[Ru(CN)6]*- 


[Fe(CN)6]^- 


11 


13 


1.2 


sihca bead 


APTES 


[20] 


La3+ 


[Ru(NH3)]3+ 


560 


730 


1.3 


sihca bead 


poly-L-lysine 


[16] 


[CoCl2H30N8]'+ 


La3+ 


190 


120 


1.6 


sihca bead 


poly-L-lysine 


m 


[CoCl2H30N8]^+ 


La3+ 


170 


180 


1.1 


sihca bead 


chlorosilane 


m 


La3+ 


[Ru(NH3)]^+ 


130 


210 


1.6 


sihca bead 


chlorosilane 


[17j 


[CoCl2H30N8]'+ 


[Ru(NH3)]^+ 


210 


450 


2.1 


poly-L-lysine 


silica bead 


[20] 


[Ru(CN)6]^" 




22 







charge inversion is a general equilibrium effect and that cq 
depends very sensitively on Z but lacks dependence on 
the chemical structure of the ions agree with our mea- 
surements. 

A quantitative test of the theory is possible. Equa- 
tion ^ has two unknowns, (Tbarc and A/i'^. From the con- 
secutive measurements on [Fe(CN)6]'*" and [Fe(CN)6]'^~ 
in Fig 121 we extract values of (Tbarc = -1-0.45 e/nm^ and 
A/i° == 1.4fcr, assuming that A/i'^ is the same for both 
charge states of the ion. The corresponding values of 
^ic are 9.4 and 5.8kT for Z = 4 and 3, respectively. 
This indicates that specific interactions are negligible and 
that ion correlations are the dominant mechanism behind 
charge inversion in this system. 

The same calculation for the APTES measurements in 
TablelHyield values of CTbarc = +0.2 e/nm^, Ajj." = S.O/cT, 
and fie — 5.8 and 3.5kT ior Z — 4 and 3, respectively. 
This suggests that in this case specific adsorption plays a 
larger role. The difference between the two surfaces may 
occur because the value of fic for APTES and Z = 3 ions 
corresponds to the lower end of the range of validity of 
Eq. |j2l. In addition, the surface charge was modelled as 
being uniformly distributed, whereas real surfaces con- 
sist of discrete chemical groups; the relative importance 
of this disorder should be greater for APTES with its 
smaller value of (Tbarc- 

Taking A^° = lAkT and cq = 200 fiM 
for [CoCi2H3oN8]'^+ screening silica gives ffbarc — 
—0.4 e/nm^, in agreement with commonly accepted val- 
ues 113]. 

These experiments are among the first systematic steps 
toward understanding the fundamentals of screening of 
real surfaces by multivalent ions. Specific binding does 
not provide an adequate explanation for our observations. 
An alternative description based on ion correlations pro- 
vides qualitative and semi-quantitative agreement with 
observations. In the future, measurements using electro- 
static gating will allow tuning the surface charge density, 
permitting further quantitative tests of the theoretical 
predictions. 
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